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The quantum evolution equations for the field expectation value are analytically solved to cubic 
order in the field amplitude and to one-loop level in the model. We adapt and use the renor¬ 
malization group (RG) method for such non-linear and non-local equations. The time dependence 
of the field expectation value is explicitly derived integrating the RG equations. It is shown that the 
field amplitude for late times approaches a finite limit as 0{t~^^^). This limiting value is expressed 
as a function of the initial field amplitude. 


I. INTRODUCTION AND MOTIVATIONS 


An important research activity develops presently on the non-equilibrium dynamics of quantum field theory with 
multiple physical motivations 10-^,|j. They arise in the reheating of inflationary universes, the eventual formation of 
disordered chiral condensates , the understanding the hadronization stage of the quark-gluon plasma as well as 
try i ngi to understand out of equilibrium particle production in strong electromagnetic fields and in heavy ion collisions 

0 , 0 . 

A common feature in all such phenomena is the presence of quantum fields with large amplitudes. That is, they 
imply the quantum field dynamics when the energy density is high. That is to say of order m~'^ times a large number, 
where m is the typical mass scale in the theory. [For example, where A is a coupling constant]. The usual 

S-matrix calculations apply in the opposite limit of low energy density and since they only provide information on in 
—> out matrix elements, are unsuitable for calculations of expectation values. 

Our programme on non-equilibrium dynamics of quantum field theory, started in 1992 j^], is naturally poised 
to provide a framework to study these problems. Thorough and accurate numerical calculations involving large-V 
and Hartree nonperturbative methods gave a clear picture of the physical processes involved both, qualitatively and 
quantitatively |2|jj,p. Very briefly, part of the energy initially stored in the zero mode is dissipated through quantum 
particle production. Since the initial field amplitude is typically of order l/vT, this is a nonperturbative physical 
process. Detailed analytical calculations yield precise results during the so-called preheating time. Namely, before 
non-linearities shut-off the process of particle production linked to parametric or spinodal instabilities [^. 

We consider here a scalar field theory in four space-time dimensions (i.e. a standard inflationary model). We 
give the name of order parameter to the quantum expectation value of the field 4>{t). We consider traslationally 
invariant quantum states with constant (and finite) energy per unit volume. A nonlinear and nonlocal equation of 
motion determines the time evolution of To one-loop level and to cubic order in the equation of 

motion can be written as 


6 


-4>{t) / dt' ^{t') 


'to 


d^k cos[2u!k(t — t')] 

(^' 


2 ^^ 


= 0 . 


( 1 . 1 ) 


where g and w? stand for the renormalized coupling constant and mass and ujk = J- Eq.(l.l) is to be 

supplemented by the initial data (/'(to) and 4>{to)- _ 

The purpose of this paper is to solve analytically the equation of motion 0) within the renormalization group 
approach. We obtain in this way analytic expressions for the order parameter at all times. 
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Our final results can be summarized as follows. The order parameter takes the form 


(^(t) 



ri{mt) 


where 


r]{t) = R{t) cn(t\/l + R{tY,k{t)) , 

where cn stands for the Jacobi elliptic cosinus and R(t) and k{t) = ^ 2 t^+R(tY) slowly with time. [Slow compared 

with the variation of the elliptic function whose period is of order one]. 

The time dependence of R{t) is determined by the renormalization group equations. We give in eq.([^) the explicit 
time dependence. 

We find for late times, 

A-y/^ Rl^ ricosnt cos(2t + 7r/4) + sinOt sin(2t + 7r/4) - 2 \ 

R[t) — Roc — 22 I _ 3^2 TrTTWTSTr^! Tm72 ^ ) 


(n + 2)2(f2-2) t3/2 


That is, the field amplitude for late times approaches a finite limit as 0{t This limiting value can be expressed 
as a function of the initial amplitude Rq by 


F{Roc) - F{Ro) =-^ 

Lot 




Arg cosh wi — 2 


where the trascendental function F{R) is defined by eqs. (El) and ( |7.8| ). 
For small initial amplitude this reduces to 

+0(Ai?4) 


In summary the order parameter oscillates as the classical cnoidal solution with slowly time dependent amplitude 
and frequency. The amplitude decreases with time due to the energy dissipated in produced particles. These results 
are in agreement with the numerical calculations |Q,^. 


II. A PEDAGOGICAL EXAMPLE 

The renormalization group method (RG) |pl , |il , |il| improves the knowledge of the asymptotic behaviour of solutions 
of non-linear differential equations. The long time power law behaviour for non-linear equations can be found by RG 
methods [ pRljl^ ]. The RG not only resums the secular terms that arise in an usual (naive) perturbative expansion 
but also yields an uniform expansion of the solution. More, it is possible to choose between the different resumed 
approximations. The different choices are, of course, all equivalent at a given order. We will show that the choice of 
the final form is given from the choice of the counterterms that are needed in the method. Moreover, after one has 
choosen the final form of what can be called the envelope of the solution [Q , one still has a freedom in the choice of 
the counterterms and can introduce an arbitrary function expressing a kind of gauge freedom. This gauge function 
may be used to give a simple expression to the slowly varying pieces. Of course, the final solution is independent of 
this gauge function. 

In order to introduce these ideas we consider the exactly solvable example of the anharmonic oscillator with given 
initial conditions: 

y + y + ey^ = 0 
y{to) = Rq, y{to) = 0 

Where the exact solution is given in terms of elliptic functions 11 
y{t) = Rq cn - to)yll + eR^, k{^/eRo)^ , 


k{x) = 




( 2 . 1 ) 
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Starting from the zeroth order solution (for e = 0) 


yo{t) = Ro cos{t - to) 

one can easily get the first order solution by linear perturbation (or by expanding eq. & to order e) and see that it 
contains a secular term (in t — to)' 


y(t) = Ro cos(t + 0o) - e to) i?o sin(t + ©o) 

O 

R^n 


32 


cos 3(t + ©o) — cos(t + ©o) 


+ 0(e^) 


©(^o) — — 'to 

Now following the method described in ref. we introduce the renormalization of Ro and ©q: 


( 2 . 2 ) 


i?o(to) = R{T)Zl{to,T) 


©o(to) = ©(t) + Z2{to,T) 


Where r is a constant which we shall choose later, we expand Zi and Z 2 in power of e, after we we have noticed that 
Zi is a multiplicative renormalization and that Z 2 is an additive one: 


^1 

^2 


Inserting these expansions in (2.2) we get: 


1 + eoi + O(e^) 
€bi+0{€^) 


y{t) = 


Rcos{t + Q) + e Rai cos{t + Q) — e R bi sin(t + ©) 

3 3 

— e -{t — t) R^ sin(t + ©) — e -(r — to) sin(t + ©) 


+ e 


i?3 


cos3(t + ©) 


cos(t + ©) 


0(e2) 


(2.3) 


Now we want to remove the term that grows with to, because for a fixed (eventually large) t, this is a secular term. 
Although there are no ultraviolet neither infrared divergences, which are typical of quantum field theory, subtracting 
such terms is the analogue of eliminating the divergent terms in a quantum field theory calculation (to <-*■ In A), where 
A is the ultraviolet cutoff. We call minimal subtraction scheme the one with the following counterterms 


a"® = 0 
tr = I (fo - r) 


Requiring that y{t) do not depend on r (as bare correlation functions do not depend on the renormalization scale in 
quantum field theory), we obtain: 

0 = ^(t) = R(t) cos(t + ©(t)) + eR 
dr 

Using the fact that for a fixed r, cos(t + ©(t)) and sin(t + ©(r)) are independent functions; one can see that they are 
both solution of the zeroth order equation. This yields the renormalization group equations, (which can be obtained 
using different but equivalent means 0@): 

R = 0 + C>(e2) 

e = elR‘^ + Oie‘^) (2.4) 

8 


|r(t)2 - 0 


sin t + ©(T) +0(e") 


Now we consider that R and © are solutions of the RG equations, hence y(t) is r-independent and we can choose r 
arbitrarily. We choose t = t and we recover the same result as expanding the exact elliptic solution (see ref. ffl). 
(This approximation can be derived by other methods, see ref. ||l4||): 
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y{t) = Rq cosu}{t — to) + € 
w = 1 + e -i?o + C>(e^) 


32 


cos3u}{t — to) — cosuj{t — to) 


0{e^) 


(2.5) 


Here the secular term have been resumed [compare with eq.(2.2)]. 
We could have chosen other counter-terms, for instance: 


,MS 

R2 

h 


MS 


1 

^ 32 


— 1 -I- cos2(t -I- 0) — g{T -I- 0) sin(T -|- 0) 
sin2(T -I- 0) — g(r -I- 0) cos(t -|- 0) 


where g{u) is an arbitrary function. One has to not ice that the terms in sin2it and cos2u are introduced to remove 
the term in cos3(t -I- 0o) in the hrst order solution (^^). 

When we insert these values in (2.3) we get: 

3 

y(t) = R cos{t + 0) + e - (r — t)R^ sin(t -|- 0) 

8 

i?3 


32 


cos3(t -I- 0) — cos(t -I- 0) cos2(r -|- 0) -I- sin(t -|- 0) sin2(r -|- 0)-|- 


g{T + 0)(sin(T -1- 0) cos(t -f 0) — cos(t -|- 0) sin(t -|- 0) 


0{e^) 


The RG equations take now the form: 


R = -e 


32 


2sin2(r -|- 0) -|- 


R 


dr 


gir -I- 0) sin(T + 0) 


0[e^) 


Q = e-R +e—{2cos2{T + Q)- — 


32 


g{T + 0) cos(r -f 0) 


Oie^) 


If R{t) and Q{t) are solutions of these last RG equations, y{t) is r-independent, but might depend on g{u). In fact 
this is not the case, if one puts t = t and expands the final solution to see how it depends on g{u) one notices that 
at this order all the terms containing g{u) cancel. In other words, the hnal solution do not depend on the choice of 
this function g{u). The arbitrariness of the function g(u) is like a gauge freedom in the present context. We can use 
this gauge freedom to simplify the RG equations. 

In the present case a good choice is: 


g{u) = —2 sinu 

With this value of the gauge function the RG equations become: 

R = 0 + 0{€^) 

0 = € ^R^ + € ^ cos 2(t -I- 0) -I- 0(€^) 

8 8 

It is clear that if we again choose t = t, all terms will vanish except the first one: 

y{t) = R{t) cos(t + 0{t)) + 0{e^) 


We finally get: 


y{t) = Ro cos 
3 


j^2 

uj{t — to) + e— sin 2uj(t — to) 
lb 


O(e^) 


oj — I -|- T C^(c ) 


Expanding this last result in powers of e we recover the previous equation (2.5). 
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III. OUT OF EQUILIBRIUM QUANTUM FIELD THEORY 


Let us consider the self-coupled scalar field theory with Lagrangian density 

£ = _ ^$4 

The time evolution of the order parameter defined as the quantum expectation value, 


(3.1) 




(3.2) 


is governed by non-linear and non-local equations. We consider translational invariant quantum states. 

As shown in ref. [^, to order g (one-loop level) and keeping the first non-linear contribution in the field amplitude 
4>{t), the equation of motion for takes the form: 

^{t) + dt' ^{t') J —— = 0 • ( 3 - 3 ) 

Here g and stand for the renormalized coupling constant and mass, respectively. We have chosen a finite time to 
as the beginning of the self-coupling interaction as in ref. 

It must be noticed that the last term in eq. (^.3|), (the ‘dissipative’ contribution), has a non-Markovian (i.e. 
memory-retaining) kernel. Secondly, the equation is non-linear the field amplitude (/'(t). Moreover, as shown in [^, 
a.V(j) term in the equation of motion cannot reproduce the dynamics. 

Now, we introduce dimensionless variables and replace thy tjm 


and eq.(^) becomes 





i) + g + g^ + Xh[rf\ = 0 , 


(3.4) 


where 


h[g{t)\ 

Kit) 


Vit) [ vit')vit)Kit' -t)dt' 

Jto 


1'°° k^dk cos 2ujt 

Jo 27r 


and 


1 + k^ 


(3.5) 


Since to is the beginning of the out of equilibrium evolution the initial conditions for t 


to are given as in ref. 


r]{to) = Ro and 77 (^ 0 ) = 0 


In order to understand the late time behaviour of the order parameter, we shall solve this non-linear equation 
Using the renormalization groups methods | p^ , pU 181 we will derive this behaviour in term of the initial conditions, 
up to the first order in A and third order in the amplitude. 


IV. THE NAiVE PERTURBATIVE EXPANSION 

We start by seeking a standard perturbative solution in powers of A, 

r]{t) = r]o{t) -hXm(t) -\-0{X^) , 

with goito) = Ro and 770 (^ 0 ) = 0 . (4.1) 

We have then to solve to zeroth order: 


??o + ?7o + ??o = 0 , 
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r]o{to) = Ro and 770(^0) = 0 . 


(4.2) 


The first order correction obeys, 

m +m+ + h[r]o] = 0 , 


? 7 i(^o )=0 and 771 (^ 0 ) = 0 . 


(4.3) 


We now introduce the function: 


with 


fit, R)=R cnitVl + R'^. HR)) 

HR) = , ^ ■ 

x/2(l + i?2) 

It is straightforward to show [|3| that /(t, R) obeys eq.(4^) for to = 0. 

f + f + f=0, 

fiO)=R , /(0) = 0. 

So we have to zeroth-order the (purely classical) solution: 

Vo{t) = f{t-to,Ro) ■ 


(4.4) 


(4.5) 


In order to find the one-loop quantum correction rji, we have to solve the linear differential equation (4.3). We define: 


f,it,R) = ^{t,R) 


f2it,R) = ^it,R) . 


As fit) is a solution of (|4.2D and /i(t) and / 2 (t) are both derivatives of fit), /i(t) and / 2 (t) are solutions of the 
homogeneous part of ([4.3|). Using the Green’s function method, we can write the solution of eq.(4.3) as 


771 (t) = - 

' ^ ^ W'o 


f 2 it-to) f flit'- to)hit')dt'- flit - to) f f 2 it'- to)hit')dt' 
J to J to 


(4.6) 


Where hit) stands for ti[77o(t)] and where Wq is the Wronskian of /i(t, R) and / 2 (t, R)'. 

Wo = W[fiit, R), / 2 (t, R)] = -Roil + Rl) ■ 

We now have the first two terms of the p ertu rbative expansion of rjit) in powers of A . Unfortunately, they contain 
secular terms that make the expansion (© useless for t > A To see this, one can Fourier expand the periodic 
functions fit) and /i(t) (see [O) as: 


fit, R) = ^fn cos( 2 n - l)nt 

n—0 

00 

fiit,R) = -n E i2n — 1) fn sin(27T, — l)Ut 


n—0 


n = 


-yrr^, fn = 


2ttR 


1 - 1/2 


q = e 


-■KK'(k)/K(k) 


(4.7) 


Kik) •"* k Kik) > 

One can notice that the function / 2 (t) is not periodic in t. However, the function foit) defined as 

fsit) = f2it)-t ^fiit) , 

is a periodic function of t with the Fourier expansion 

00 

Hit) = E cos( 2 n - l)nt . 

n—0 

Here ' stands for 

This shows that Hit) contains a secular term, that will contribute through eq.(4.6), to the secular terms in 771 (t). 
The integrals of /i(t — t')hit') and / 2 (t — t')hit') in eq.(4.6) yield secular terms. The second expression produces a 
second order secular term (growing as t^). 
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V. RENORMALIZATION GROUP EQUATIONS 


In order to analyze the late time behaviour of the parameter, we want to recast ri{t) in the form: 

Vit) = f{t + Q{t),R{t)) , 

where the t dependence of 0 and R is slow. That is, 0 = 0(X) , R = 0{X). Since f{t) is an amplitude R multiplied 
by an oscillatory function (the elliptic cosine), R{t) will be the envelope of the solution. 

In order to do that, we write the parameters Rq and 0o(io) = ~to in terms of new parameters 0 and R through 
(finite) renormalizations: 

Ro{to) = R{t) Zi{to,T) 


0o(to) = Q{t) + Z2{to,T) 

Here r is a (new) constant parameter and Z\ and Z^ have the following expansion in power of A, which ensure that 
both R and 0 are of order 0{X): 

= l + Aai(to,T) + 0(A2) 


Z 2 — 0 + A &i [to, t) (D[X ^) 

The requirement that r]{t) = r]o{t) + Xr]i{t) + 0{X'^) takes the form f{t + 0(t), R(t)) up to 0(A^) fixes the (finite) 
counterterms oi and bi up to an arbitrary ‘gauge’ function g(r). We find using eqs.([4.5|) and (|4.6|), 


ai{to,T) = ^/i('^ + 0) ~ ^ fi{t' + e)h{t')dt' 
bi{to,T) = + 0) - ai(to,T) (r- to)^ 


1 

W 


(r-to) / flit'+ Q)h{t')dt'- f f2it'+ Q)hit')dt' 

Jtn Jtn 


'to to 

Where we have written 0 for 0(r), R for R{t) and fi{t) for fi(t, Rir)). We find, 

Tjit) = fit + e,R)-X^ [flit + e)MT + 0) - hit + 0)/i(r + 0)] 

+A + 0)/i(i + 


+ 


w 


hit + 0) + (t — to)—flit + 0) 


flit' + e)hit')dt' 


-fiit + e>)J^ hit'+ Q)hit')dt'^ +oix‘^). 

Since ry(t) is clearly r independent: 

drj 


dr 


W = o , 


vt 


Imposing eq.(^) to eq.(^) yields the renormalization group equations to order A. 


= ^r\Ki)hit + ^i't))- 


IT 

A 


dt 


IT 


m = -^AHt)f2it + eit))-- 


dt 


git) flit+ Qit)) 
git) hit + Bit)) 


(5.1) 


(5.2) 


(5.3) 


In the calculations, we used the fact that 0 and R are of order 0(A) and the line ar independence of fi and h- 
Since r is arbitrary and 77 (f) is t independent we can set r equal to t in eq.(5.1). We thus find, 

77 (f) =/(f + 0 (f), i?(f)) + 0 (A 2 ) . 

We have thus succeeded to resum the secular terms. 
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VI. FIXING THE ‘GAUGE’ FUNCTION 


We introduced an arbitrary function g^r) in the counter terms oi and bi. This function appears in the renormal¬ 
ization group equations. The solution r](t) must be independent of the choice of gir). This is easy to prove to this 
order. One can choose the freedom in g to give to 0 or i? a specific form. It is clear that 0 is related to the finite 
mass renormalization and so has to be bounded for large t. 

We now study 0 from eq.(5.3). We have to compute 

/i(t) =/(r-b 0) J dt' ff{t' + Q) J dfi{uj) cos — t) . 

Where the integration measure for uj is 


dpL{iu) = dk = du 


27ra;3 


2'ku)'^ 


We introduce 


I{t) = / dt' ff{t' + Q) / dfi{ui) cos2ui{t'— t) 


such that h(r) = /(r -b 0) /(r). 

Using the expansion of the Jacobian elliptic function (see ^l|) and the notations © we have 




K{k) 


nq 


E m 

T~r 


^ l-q- 

n—l 


■ COS nQt 


( 6 . 1 ) 


and 


OO 9 n 


ffi{t -b 0) = -20^ ^ ^ sinnn(t -b 0) 


n—l 


where K{k) and E{k) stand for complete elliptic integrals of first and second kind, respectively. 
Then it is straightforward to obtain that: 


( 6 . 2 ) 


i[t) = 

^ I _q2r 

n—l 


dg{uj) 


cos nil{t -b 0) — cos(nn0 — 2ujt) 


2u; + nfl 


cos{2ujt -b nflQ) — cos nn{t + 0) 


2u} — nVL 

We study the behaviour at large t of each term of I (t). For the first integral we trivially get 


/■ cosnr2(t-b0) 

J — 2uj + nn — ^ cosnVt{t + 0) 

1 


C« = 


We find for the second integral 


- / dg{uj) 


cos(nfl0 — 2uit) t- 


2uj -b nVl 

The last integral can be written as 

cos(2cli< -b nflQ) — cos nid{t + 0) 


2TTio^ 2io + nQ 


/it cos{2t — nflQ + IT / 4) 1 

8 (2 + nfI)t3/2 


dgico) 


2io — nQ 


= c^\t) cosnn(t-b 0) — s)f^(t) sinnn(t-b0) 


Rf - Fh 


(6.3) 


(6.4) 


Where, 
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Since the final resumed solution do not depend on g{t) and since ri(t) is finite we may choose g(t), in order to find a 
simple expression for R{t), as a solution of the equation 

Idzv{t)ffi{t + 0) = ^ [9it)fi{t + 0)] 

That is, 

+ 0) dt' . 

Now the renormalization group equations become 

0(0 = “ ^div{t) ffait + 0) + Iregit) //2(< + ©) 

“ It + ^git)flit+ 0)'^ 

m = ^Iregit) fflit+ e). ( 6 . 6 ) 
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VII. THE ORDER PARAMETER BEHAVIOUR FROM THE RENORMALIZATION GROUP 

EQUATIONS 


W e have now the renormalized group equations (6.6) an d w e mus t rem ember th at the field equations of motion 
(3.4) are valid up to 0{Rq). We have, using the expansion (^) and (^) in eq.(^) 


2 Tl 


Idiv{t) - X! 1 _ „2n 
n—1 ^ 

Ireg(t) = I(t) - Idiv{t) 

We now introduce the following notations: 

1 


Cn COS nVL{t + 0) — 5'„ sin nVL{t + 0) 


— 


, F\R) = 


K{Rf 


R 


^ ^ q{^B)2 (1 + i?2)2 

With this notations the renormalization group equation ( |6.6| ) for R{t) can be written as: 

tt^A { f \ 

F'{R)R=—^ ^ 7rp(a;„) ( cos2(wm - w„)t - cos2(a;m + Wr!,)M 

n,771=1 t ' 


(7.1) 


(7.2) 


+reg 


1 


1 


U! — LOn LU FLUn 


sin 2{ojjn + uj)t — sin 2(w — ujm)t I did{uj) 


Where whe re by ‘reg’ we me an the i ntegral reg ularized by subtracting the appropriate cou nterterm contained in Cn 
[see eq.(^)]. We used eqs.(|^-6.3) and (M) and we have set to = 0 since the eq.(^) is invariant under time 
translations. Moreover, at this order of A we can set 0 = 0. Integrating this last expression with respect to t yields: 


77 /n ^ AtT® j f sin2{uJn - UJm)t Sill 2 {uj m + UJn)t 

F{R{t)) - F(Ro) =—— > anam<Trp{(^n)[ -;- 

4 I \ 

71,771 = 1 ^ ^ 

/ oo 

dp.{uj) 


UJ — UJr). UJ -\- Ur, 


1 — COS 2(u; + Um)i 1 — cos 2(u; — Um)i 


U) -|- U)r 


(7.3) 


The terms which look singular for n = m must be understood in their limiting value n ^ m, which is always regular. 

One should notice that, unlike the other terms of the integral, the last term gives (for n equal to m) a secular term, 
which is of course canceled by the last first 2'!Tp{ujn) t: 


/'°° 1 - cos2(w - w„)t 

/ -7-— 

/l (Ul-UJn)^ 


/ C 


7 ^+ 0 ( 4 ) 

[UJ - UJnr t 


The equations of motion ( |3.4[)- (p7^) are valid to order i?Q and q = 0{k‘^) = 0( Rn) fo r small Rq- Therefore, we can 
keep the terms in the sum (|7.!^) where n and m are equal to 1. We find from eq.(|7.3|) after non-trivial cancelations, 


F(i?(t)) - F{Ro) = 




-Kp(uJl) ( 2t — 


sin 4a; it 
2a; 1 


/ OO 

dp{uj) 


Ul — Uli OJ + OJl 


1 — cos2(a; + a;i)t 1 — cos 2(a; — a;i)t 


a; + a;i 


a; — a;i 


(7.4) 


where a;i = Q/2 = 2 K(k) ^^ 

In order to find the field amplitude for t —> oo, we average over an oscillation period 


foo 

FiR^) - F{Ro) = -— v.p. / duj 






(a; + a;i)2 (a;-a;i)2 


The last integral gives: 
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F{R^) - F{Ro) = 



ujf + 2 
- 1 


Argcoshwi — 2 


(7.5) 


One can easily see that the expression between the brackets is positive. Since F'{R) is positive i?oo is smaller than 
Rq. This result was expected since we have a dissipative behaviour. Namely, part of the energy initially in the zero 
mode is spent in particle production. 

More explicitly, 


F{R) = r 

J oo 


K{k{R)f RdR 

q{k(R)Y (l + i?2)2 "A/vi diky ■ 




(7.6) 


where k{R) is defined by eq.(4.7|). 

For late t, R{t) — i?oo is small and we can write 

F{Rit)) - F(i?oo) = F'iR^)iRit) - i?oo) . 

We can then evaluate F{R{t)) for late t from eq.(|7.4|) using eq.(|6.4D, 

A(i?(<)) - F^(i?oo) = TT® 


/■°° cos(2t'+ 7r/4) sinflt' , 

at 


(2 + 


We find equating both results the behaviour of R{t) for late t: 


Ay^ R^ ricosAlt cos(2t + 7r/4) + sinlAt sin(2t + 7r/4) - 2 \ 

R(t)- R^ = - - — - ^ U(t ) 


32 l-3i?2 


(0 + 2)2(f4-2) t3/2 


(7.7) 


Eqs.(^) and (^) define i?oo as a function of Rq and A. 

In the present case it is very convenient to express the elliptic functions in powers of the elliptic nome q = e~'^^ 
since q is always small in our context. Namely, 0 < g < e~'^ = 0.0432139... for 0 < i?o < oo. We have for example, 

K(k)=^- [l + 4g + 4g2 + 0(g4)] 


k'^ = 16g [1 - 8g + AAq^ + Oiq^] . 


The function F{R) becomes 


F{R) = -- 


^6 r 


-81ogq- 12g + C>(g2) 


+ C . 


(7.8) 


where the constant C is given by C = ^ [l + 87re ^ — 12e + 0(e 377)j = 5738.9... 

Eq.@ and ([7.8[) indicates us that qoo = q{Rao) and qo = q{Ro) are related as follows. 


<?oo -9o = -cA + C>(A^) . 


where, 


c = 


8Qo 


( 171+2 


LOiy^Lof - 1 


Arg cosh 071 — 2 


[1 - 8qo + 0{ql)] 


More explicitly, for i?o 0 we find: 


That is: 



Rt 

1287r 


+ 0{Rl) 




+ 0{\Rt). 
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